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Summary 

An accurate and yet simple Meshless Local Petrov-Galerkin (MLPG) formulation for analyzing 
beam problems is presented. In the formulation, simple weight functions are chosen as test functions. 
The use of these functions shows that the weak form can be integrated with conventional Gaussian 
integration. The MLPG method was evaluated by applying the formulation to a variety of patch test 
and thin beam problems. The formulation successfully reproduced exact solutions to machine accuracy 
when test functions with C 2 continuity and an appropriate order of basis functions are used. 

Introduction 


Meshless methods are increasingly being viewed as an alternative to the finite element method [1- 
3]. Recently, a meshless Galerkin formulation was presented for beam (C 1 ) problems using a 
generalized moving least squares (MLS) interpolant [2]. In this Galerkin formulation, the derivatives of 
the weight functions used in the trial and test functions higher than the second order showed 
discontinuities (scissors) at the boundaries of the supports of the local sub-domains. Conventional 
Gaussian and other integration schemes were found to be inadequate to integrate the weak form 
accurately - very accurate and elaborate integration schemes were needed. 

In this paper an accurate and simple alternative Petrov-Galerkin formulation for beam problems is 
presented. In the present Petrov-Galerkin method simple weight functions are chosen as test functions 
and conventional Gaussian integration is used. The effectiveness of the MLPG method is evaluated by 
applying the formulation to a variety of patch test problems. 


Development of the Petrov-Galerkin Formulation 


The notation of reference 2 is used in this paper for brevity and convenience in presentation. The 
MLPG equations are 


(node) ^ + ^ (bdy) ^ _ ^ (node) _ ^ (bdy) _ ^ 

where 

{d^jw 1 , 0 1 , u 2 , Q 2 , ...} 


( 1 ) 

( 2 ) 


are the fictitious nodal values of deflections u and slopes 0, and the matrices in Eq. (1) are defined as in 
Eq. 3 5 -3 6(g) of reference 2. 
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The MLPG equations are derived using a weighted residual weak form of the governing equations. 
The trial functions used for beam problems are derived using the generalized MLS interpolation [2] as 

«(*) = xf “ Vf 0 (*) + ® V/ (0) oo] ( 3 ) 
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|a]=p ' wP + pJwP x . 

In Eq. (5) P is an (n,m) matrix and w is an (n,n) matrix defined as 


[ p ]=[pOi) PO 2 ) ••• pO„)] T , 


where x = x — x t , and 


t f 2 m - 1 

p (x) = \ 9 x, x , . . . x , 


T . - dp (x) r m-2 1 

p x (x) = — = 0, 1, 2x, ... (m-l)x 

dx L 


with (tw- 1) as the order of the basis function p(x) used in the MLS approximation. Three different 
weight functions w t (x) are considered: 


w,-(X) = 


1- (x-x,.) 2 !R- T if Wx-x^R, 

lL J 


if \\x-Xi >£,• 


with a = 2, 3, and 4. Note that the extent of the influence of the trial function is controlled by distance 
R- 



Several different test functions v 2 (x) are considered: 


v<(*)= 


(X - Xj ) 


/r: 


o 


if |x — Xj | < R 0 
if llx - x,- II > R„ 


with fi = 2, 3, and 4, 
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with 4 = || x-Xi ||, and 
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if 0 < d t < R 0 


if dt>R 0 
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The derivatives of the test functions are evaluated at the center (d t fR 0 = 0) and at the end points 
(di/R 0 = 1) as 
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The test functions are then C Y continuous up to the order y where y = m j) (See reference 4, pp. 

63-64). With these definitions, the test functions in Eq. (10) with fi = 2, 3, and 4 are C 1 , C 1 , and C 3 
continuous, respectively. Similarly, the spline functions in Eqs. (11) and (12) are C 1 continuous. Note 
that the lengths R t and R 0 in Eqs. 9-12 are user defined in the MLPG method. 

As the test functions chosen (Eqs. 10-12) are not the same as the trial functions, the current 
implementation is a Petrov-Galerkin method. This is in contrast to reference 2 where a Galerkin 
method was used. The Galerkin method led to discontinuities (scissors) at the boundaries of the 
supports of the trial functions. Due to these scissors, elaborate numerical integration schemes were 
needed to integrate the weak form accurately. The domain of dependence Q s was subdivided into 
subregions depending on the ends of the support domains of various trial functions. A 10-point 
Gaussian quadrature was used in each of these subregions to integrate the weak from accurately. In the 
current implementation with the Petrov-Galerkin test functions, the integrands were examined and no 
discontinuities were found in the domain of integration. A single domain for the entire compact 
support domain (Q s ) was used. A standard 8-point Gaussian quadrature rule was found to be sufficient 
to integrate the weak form very accurately. 



Beam Configurations and Models 


A beam of constant flexural rigidity El and a length of 41 is considered. Six models with 5, 9, 17, 
33, 65, and 129 nodes uniformly distributed along the length of the beam are considered. Figure 1 
shows a typical 17-node model. The distances between the nodes (A//) in these models are 1.0, 0.5, 
0.25, 0.125, 0.0625, and 0.03125, respectively for the 5-, 9-, 17-, 33-, 65-, and 129-node models. Three 
types of basis functions, quadratic basis (1, x, x 2 ), cubic basis (1, x, x , x 5 ), and quartic basis 
(1, x, x 2 , x } x 4 ) are used. The local coordinate approach of reference 5 is used in evaluating the shape 
functions and their derivatives. System matrices in Eq. (1) are developed with these parameters. 
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Figure 1: A 17-node model of the beam 

Numerical Evaluations 

The current MLPG formulation was evaluated by applying the formulation to simple patch-test 
problems. Three typical problems were considered. 

du 

1 . u(x) = c 0 , 0 = — = 0 ; Rigid body translation 

dx ( 14 ) 

2. m(x) = C!X, 6 = Ci ; Rigid body rotation 

2 

3. u(x) = c 2 x , 0 = 2c 2 x ; Constant - curvature condition 

All three of these problems satisfy the governing differential equation exactly and as such are exact 
solutions. The deflection u and the slope 0 corresponding to problems 1, 2, and 3 were prescribed as 
essential boundary conditions (EBCs) at x=0 and x=4 /. With these EBCs, the beam problem was 
solved using the MLPG method in Eq. (1). When ( R 0 /l - 2A) and (R t // = 8A) were used, each of the 
5-, 9-, 17-, 33-, 65-, and 129-node models reproduced exact solutions at all internal nodes to machine 
accuracy thus “passing” the patch tests. The analysis was repeated with each of the test functions in 
Eqs. (10-12), and the MLPG algorithm passed the patch tests for all cases. 

As mentioned previously, the parameters (R 0 / 1) and (R t /l) in the MLPG method are user- 
controlled. The previously mentioned lengths (R 0 / 1 - 2 A) and (R t / / = 8 A) were used at all nodes of 
an N-node model, except at node 2 and node N-l (see Figure 2). For these nodes (. R 0 /l - A) was used 
to ensure a symmetric Q s . Note that with these assignments of (R 0 / 1) the test functions for all interior 
nodes have symmetric Q s configurations. As shown in Figure 2, no asymmetry is introduced at nodes 1 
and N as exactly half of their test functions are used. When these symmetries were violated, the MLPG 
method failed the patch tests. 

As the models are refined, the value of (A 11) decreases and thus the size of Q s and the extent of the 
trial functions also decrease. For finer models, i.e. for the 33-, 65-, and 129-node models, when 
8 A < R t f l < 16A the MLPG method yielded very accurate results. However, when R t / 1 > 16A the 
MLPG method failed the patch tests for these models. When R { / 1 > 16 A the trial function is too 



diffused and the size of Q s (R 0 // = 2 A) is too small in comparison to (R t / V). The combination of small 
Q s size and large (R { / T) are apparently incompatible. While the finer models performed well over a 
large range of (R t / /), the coarser models performed well in a much smaller range of (R t /l). For good 
performance, (R { / 1) needed to be around 8 A but less than 98% of the total beam length. 
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Figure 2: Q s definitions for various nodes 

The MLPG method was applied to beam problems with mixed boundary conditions. The first 
problem considered was a cantilever beam with a concentrated moment at the free end (i.e. M = M 0 at x 
41). The exact solution for this problem is u = MqX 2 / 2EI and 0 = MqX / EI. For all trail functions 
considered, the MLPG algorithm reproduced the exact solution when the test functions in Eq. (10) with 
fi = 3 and 4 and when the 4-term spline function in Eq. (12) were used. In contrast, the algorithm failed 
to reproduce the exact solution when the test function in Eq. (10) with fi = 2 and the 3 -term spline 
function of Eq. (11) were used. This example suggests that test functions with at least C 1 continuity 
and with m x > 2 (see Eq. 13) are required for the MLPG algorithm for beam problems. 


The second problem considered was a cantilever beam with a tip load. Since the exact solution for 
this problem is a cubic in terms of the x-coordinate of the beam, all six models with cubic basis 
functions and a test function with C 1 continuity and with m x >2 reproduced the exact solution to 
machine accuracy. 


The third problem considered was a simply supported beam subjected to a uniformly distributed 
load. Using symmetry, half of the beam was modeled. Since the exact solution for this problem is 
quartic in terms of the x-coordinate of the beam, the MLPG method with a cubic basis function did not 
reproduce the exact solution. Error norms defined as 
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were computed at g uniformly spaced points along the beam. A value of g = 200 was used. The norms 
HEJI 2 and \\E m \\ 2 are presented in Table 1. The numbers in Table 1 were computed using a 20-point 
Gaussian in each of the single compact support domains £2 S . As expected, all models yielded accurate 



solutions (within 4% for u, 6 , and M). When the order of the basis function was increased to quartic, 
the MLPG method reproduced the exact solutions (for u, 0, M, and V) to machine accuracy. 

Table 1: Error norm ||E|| 2 for a simply supported beam subjected to a uniformly distributed load 
with cubic basis used in the MLPG method. (Trial function using Eq. (9) with a = 3 and test 


function of Eq 

1— l 

o 

£ 

S' 

II 

4- 

Error norm 

Number of nodes in the model | 

5* 

9 * 

17 

33 

65 

129 


0.1662e-l 

0.1306e-2 

0.4573e-2 

0.3829e-l 

0.1742e-l 

0.2368e-l 


0.2774e+0 

0.1057e-l 

0.1704e-l 

0.3680e-l 

0.1763e-l 

0.2340e-l 


* Ri /l = 3.9 


The above problem demonstrates the following phenomena. When the order of the basis function 
equals the order of the exact solution, the previously stated 8 -point Gaussian in a single Q s is sufficient 
to integrate the weak form very accurately. When the order of the basis function is less than the order 
of the exact solution, a higher order integration rule (such as a 20-point Gaussian) is needed to obtain 
accurate results. Note that a 20-point Gaussian in a domain Q s is preferable over the sub-domain 
integration required in the Galerkin formulation. 

Concluding Remarks 

This paper presents an accurate and yet a simple Petrov-Galerkin formulation for solving beam 
problems. In the formulation, simple weight functions are chosen as test functions. These test functions 
must have at least C 1 continuity and zero first and second derivatives at the end points for the algorithm 
to pass the patch tests. When appropriate test functions are used, the weak form can be integrated with 
conventional Gaussian integration. The MLPG method was evaluated by applying the formulation to a 
variety of patch tests. Weight functions used to develop the trial functions need not have C 1 continuity; 
C° continuity may be sufficient. Additionally, the weak form requires at least a quadratic basis, and a 
cubic basis function is preferable. The formulation successfully reproduced exact solutions to machine 
accuracy for thin beam problems when an appropriate order of basis function was used for all 
combinations of trial and C 1 continuous test functions considered. 
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